Theory of spatial mode competition in a fiber amplifier 
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Theory of monochromatic wave field amplification in a waveguide array based on expansion of 
the wave field in terms of guided array modes is developed. Equations for the expansion coefficients 
include cross-modal gain, which completely changes behavior of the amplified wave field. Analysis 
of two-mode amplification reveals new unusual features in characteristics of two-mode amplifier. 
Instead of unlimited growth of both modes for incoherent fields, one of the modes grows with 
no limit and suppresses the lower-power mode. Effects associated with the cross-modal gain are 
illustrated analytically on a system of two thin parallel planar waveguides. Conditions are found 
when a mode with lower gain can become the dominant one at the output of amplifier. 

PACS numbers: 42.65.Wi, 42.55.Wd 



I. INTRODUCTION 

Wave field propagation over an array of parallel waveg- 
uides is up to now a subject of numerous theoretical and 
experimental studies, mostly for optoelectronic devices. 
In particular, non-linear refraction of materials can lead 
to a wide variety of non-trivial physical phenomena in 
the waveguide arrays, thus making them very appealing 
for research. Lattice solitons [Ij and non-linear optical 
switches Q are among such phenomena. An array of 
nonlinear waveguides in a multi-core fiber (MCF) laser 
Q is another example very important for the laser de- 
sign. The MCF laser construction is an array of single- 
mode cores doped by rare-earth ions in a common pump 
cladding. Such a design allows one to enlarge a combined 
mode area, thus diminishing the fiber length. Single- 
mode cores can be arranged also in a form of active de- 
fects in photonic crystal fiber lasers [i| . The MCF design 
has an advantage over single-core large mode area con- 
structions only if the fields in the cores are phase locked. 
Phase locking can be provided by an external optics or 
by diffractive exchange of fields between the cores. 

The purpose of this article is to analyze the mode com- 
petition in the multi-core fiber amplifier. This problem 
seems to be important far outside the fiber laser the- 
ory because it can be considered as a general dynamic 
problem for the Schroedinger equation with a complex- 
valued nonlinear potential. The nearest quantum me- 
chanical analogy for the optical wave field propagation is 
wave function temporal variation in an array of potential 
wells. 

The method commonly used in describing the guided 
field evolution in passive systems is a coupled mode the- 
ory Q (CMT), in which the wave field is presented as 
a sum of sole waveguides modes with unknown coef- 
ficients being a function of the propagation distance. 
These coefficients satisfy equations, which are derived by 
perturbation method. CMT is intuitively evident and 
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can be easily extended to new physical effects jf|. The 
CMT for a pair of nonlinear waveguides was analyzed in 
This theory was generalized by Hardy and Streifer 
who introduced non-identical waveguides and vector 
field effects. Phase self-synchronization at high intensity 
1 GW/cm 2 j and interaction of wave fields with dif- 
ferent polarizations and frequencies in systems of 2 and 
3 evanescently coupled waveguides were topics of earlier 
work, which was summarized in Q. 

Recently the CMT approach was implemented for 
analysis of dynamic stability for MCF configurations with 
a separate multi-channel coupler. This work was inspired 
by experiments [11]. Using model coupling matrix, the 
authors [l(| have explained the robustness of the phase 
synchronization in such systems provided special asym- 
metry is introduced into the construction. 

An MCF laser with 7-core hexagonal structure has ex- 
hibited experimentally [l2j the far field patterns typical 
for phase-locked operation at power level more than 100 
W. An analogous 19-core fiber amplifier [l3| has achieved 
20 dB gain at near to diffraction limit beam quality. The 
calculation based on the CMT predicts [l4| that the ob- 
served self-synchronization stems from intensity depen- 
dence of the refractive index of rare-earth doped fiber. 
However, direct numerical modeling based on 3D beam 
propagation method [l5[ (BPM) has clearly shown 
that the gain non-uniformity in the system is the major 
factor responsible for the effect, while the role of refrac- 
tive index nonlincarity is subsidiary. 

An important feature of the MCF amplifier is that the 
field distribution in a transverse plane is strictly deter- 
mined by the index profile, while distortions induced by 
gain being negligible. This fact approves an approach 
based on wave field expansion over passive structure 
modes, which can be found easily by a standard solver. 
In particular, such modes for the 7-core fiber design were 
found in [l7T ] by a finite element solver. 

Comparison between results of 3D BPM calculations 
(l6| and of modal analysis revealed [l8| a seeming con- 
tradiction: an in-phase mode according to mode solver 
possesses the lowest gain, while according to 3D BPM 
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approach the wave field converges to this mode, if the 
differences in phases of launched beams are within a few 
tenth of radian. To resolve this contradiction, the theory 
of light amplification in the system of parallel waveguides 
with saturable gain should be reconsidered. 

Mode competition for gain has been studied since early 
times of laser researches. Differing spatial profiles of 
modal intensities result in gain spatial hole burning ef- 
fect, which, in turn, induces instability of single mode 
lasing [l^, [2(| . Earlier studies concentrated on the com- 
petition of optical modes with different wavelengths. An 
analysis of oscillation dynamics in a two-mode laser cav- 
ity was performed in [2l|, ■ It is traditionally believed 
that this dynamics is described satisfactorily taking into 
account an effect of the modal gain saturation by the 
intensity of another mode (so-called cross-saturation, it 
was studied mainly for counter propagating modes of ring 
lasers Q). 

It is shown below that the major factor governing com- 
petition of the modes in the multi-core amplifier is a 
cross-gain (24|, which is defined as the product of the 
mode field profiles with the gain distribution integrated 
over fiber cross section area. General equations governing 
wave field propagation over the MCF amplifier are formu- 
lated in section HU The mode suppression mechanism in 
two coupled waveguides is discussed in section UIT1 A sys- 
tem of 2 delta-function type waveguides is considered in 
sections ITVl where optical modes are found, and[V] where 
simultaneous amplification of two modes in this system 
is analyzed. Amplifier characteristics are compared for 
cases of two incoherent and two coherent modes. A weak 
mode suppression mechanism by a strong one is identified 
for two-delta-function-type waveguides. The conclusion 
is made that this mechanism associated with cross gain, 
which remains unnoticed till the last time, is expected to 
exist in any system of active waveguides. 



II. BASIC EQUATIONS 

The built-in index profile of a multicore fiber amplifier 
is step-wise function with higher index regions, called 
cores, each guide single mode. To approach the goal 
to have the large mode area, small values (~ 10~ 3 ) of 
the core-cladding index difference An are used in ex- 
periments. Besides, we do not consider birefrigence ef- 
fects. It is then possible (2f| to use the scalar approx- 
imation instead of soving Maxwell's equations for the 
elecromagnetic field in the fibre. The wavefield can be 
characterized in this approximation by a scalar func- 
tion ip(x, y) exp(if3z—iut), z is the propagation distance, 
lu = kc is the oscillation frequency, k is the vacuum 
wave number, c is the speed of light, f3 is the propa- 
gation constant. The function ip could be any transverse 
component of electrical or magnetic field, obeying 2D 
Helmholtz equation, which formally coincides with sta- 



tionary Schroedinger equation in quantum mechanics: 

H^ = -/^, HV, = -(g + g)+^, (1) 

where H is Hamiltonian, f3 2 plays the role of energy, x 
and y are the transverse spatial variables, and the po- 
tential function V = — fc 2 n 2 (x,y), n is the refractive in- 
dex. Equation fl| is solved with corresponding condi- 
tions at the core-cladding boundary and at the infinity. 
This problem has a continuous spectrum of the radiating 
modes ipQ and a discrete spectrum of the guided modes 
ipj, so any field injected into the MCF can be presented 
as the sum over the spectrum of ipj plus the integral over 
the spectrum of ipQ with some coefficients .2^ . The field 
amplitudes of guided modes ipj (%> y) are real- valued func- 
tions for no-loss fiber. Besides, the guided and radiating 
modes are mutually orthogonal between each other. 

The wave field ^(x, y, z) in active fiber with gain 
g (x, y, satisfies the 3D Helmholtz equation: 

a 2 * a 2 * , j2 , , w „ 

a? + W + w + {k n - lkno9) * = °' (2) 

tiq is the cladding index. Thus, in contrast to quan- 
tum mechanics the potential in laser optics is a complex- 
valued function V + ikngg (x, y, . 

Being interested in behavior of guided modes in the 
amplifier, it is convenient to take the modes of equation 
(HJ) as a basis for expansion of the wave field: 

*(x,2,,z) = ^ Cj (z)^-e^+ J c Q (z)iP Q e^ z dQ. (3) 

3 

The non-uniform gain leads to an interaction between 
guided and radiating (leaky) modes, thus producing ad- 
ditional losses for the guided modes due to carrying away 
the energy by the radiating modes. However, transfor- 
mation of guided modes into leaky ones is a negligible 
effect for typical parameters of fiber amplifiers, as it was 
confirmed by 3D BPM calculations [26| . From the other 
side, leaky modes see the gain much lower than guided 
modes because the cladding area is much larger than the 
area of active cores. For these reasons we can substi- 
tuting expansion ([3]) into ([2]) and neglect the radiating 
modes. The gain in the cores is usually of the order of 
10 _1 cm -1 , while a typical spectral separation between 
the propagation constants of guided modes is ~ 10 cm -1 . 
Thus, variation of the guided modes amplitudes Cj due to 
amplification is slow process and it is possible to use ap- 
proximation, in which the terms d 2 Cj/dz 2 are neglected. 
Then the following system of equations for the expansion 
coefficients Cj can be derived by multiplying the resulted 
equation to ipj exp(— i(3jz) and integrating over the fiber 
aperture: 

^^ffii^EW *-^ (4) 
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where summation is made over N guided modes. The op- 
tical modes introduced are normalized for convenience as 
JJ i>j dxdy = 1. The difference of propagation constants 
is small parameter in comparison with kriQ, so the factors 
fa/ (kno) in ((4]) have been replaced by one. The applica- 
bility of this approximation was verified by 3D modelling 
of the MCF aplifier. Matrix elements gji describe an in- 
teraction between the wave field and the gain: 



9jl 



9{x,y)ipjipidxdy, 



(5) 



where integration is made over the fiber aperture. The 
coefficients gjj are the modal gains of j-th mode, while 
the coefficients gji in the sum in ^ describe the cross- 
modal gain effect. 

System (jj) of complex equations can serve as the ba- 
sis for analyzing the mode competition in a wide variety 
of MCF amplifiers. To give an idea about new features, 
which possesses the theory based on equation system (j4]) , 
let us consider the simplest situation when a waveguide 
array supports only two guided modes. In this case, sys- 
tem Q can be reduced to 3 ordinary equations for real- 
valued functions 



dPi 

dz 
dP 2 
dz 



gi\P\ + \JP1P2g12 cos <f 

giiPi + \JP~xP2gx1 cos q 

d<t> , R P1+P2 . 

T = 6/3 7WW9 12 sm ' 

dz VP1P2 



(6a) 
(6b) 
(6c) 



Here the modal powers Pj, j — 1, 2 are introduced by an 
expression Cj = y/Pj cxp(i(j>j — ifaz), <(> = fa — 4>i 1S the 
phase difference between the modal fields, 8(3 = fa — fa- 
Pj is a fraction of the wave field power carried by the 



j-th mode. Equation (f6"c| describes the behavior of the 
modal phase difference. Since gain coefficients are much 
smaller than 8(3 value, this phase difference grows along 
the propagation axis practically linearly, d4>/dz « 8(3. 
In this approximation, the problem is reduced to solving 
equations (frJaj) and (|6b|) , only. 

Spatial beatings between two modes lead to oscilla- 
tory modulation of the total field intensity along axis 
with a characteristic spatial frequency of a few 10 cm" 1 . 
This modulation induces through the gain saturation 
similar oscillations in the material gain. The resulting 
non-uniform gain profiles in different cores correlate with 
modal powers. In this case we have a tangled interaction 
between the wave field and transverse and longitudinal 
gain non-uniformities. Our purpose is to make this inter- 
action clearer using simplest waveguide configuration. 

It should be noted, that dynamic equations for the two- 
mode laser 2l| have this sort of mathematical structure 
as eqs. However, the physical content of the problem 
studied is essentially different: the authors [2l[ analyzed 
stability of two-mode regimes in a two-wavelength laser, 
while we study non-linear effects in an amplifier with 
continuous wave (cw) input signal. 



III. TWO MODE COMPETITION IN THE MCF 

To complete the system of equations (|6]) it is neces- 
sary to specify how the material gain coefficient depends 
on the total intensity. Generally, kinetics of population 
inversion varies from one type of laser to another. The 
common feature for all cw systems is so-called gain sat- 
uration, i.e. reduction of inversion and gain induced by 
stimulated emission. This effect can be qualitatively well 
described by the simplest formula g = go/ (1 + I/Isat) , 
where go(x, y) is the small signal gain, / is the wave field 
intensity, 7s at is the saturation intensity. We will nor- 
malize the intensity to the saturation intensity value. 

Even with such a simple model for gain saturation and 
taking the mode phase difference as <fi ~ 8(3 -z, the system 
is still rather difficult for an analysis. The point is 
that gji terms defined by are complicated functions of 
Pi. 2, which cannot be found explicitly for realistic waveg- 
uide structures. Nevertheless, some general properties of 
solutions to ^ can be identified. 

An important distinction of these equations is the iden- 
tity of the second terms in right-hand side of equations 
(f6"a[) and (|6b[) . For this reason it can be deduced from 
eqs. (|rja ]) - (j6b[) for the evolution of the mode power ratio: 



d (P 2 \ P2 



gi2 cos ( 



Pi- Pi 



(J) 



If (712 cos (f> is negative, then the last term in bracket of 
the right-hand side of (J7J) supports the trend the power 
ratio to grow when this ratio is greater than 1. 

The dimensionless total wave field intensity can be ex- 
pressed in a form I = P\^>\ + ft^f + 2y/PTP2'>pitp2 cos <j>. 
The cross gain variation as a function of z can be under- 
stood from an expression 



312 COS0 — 



(70^1^2 cos 4> dxdy 



1 + Pilpl + P 2 1pl + 2 Videos. 



which is obtained from ([5]). The mode amplitude prod- 
uct necessarily changes sign within the cores to provide 
orthogonality of modes. It is seen that the integrand 
generally takes on negative value with larger absolute 
magnitude when i^i^ cos <f> < 0. Therefore the quan- 
tity gi2 cos 4> is preferably negative. It can be rigor- 
ously proved that 312 cos <fi is non-positive for a system 
of two parallel waveguides possessing mirror symmetry 
An(x, y) = An(—x, y) and supporting two guided modes. 
One of modes is symmetric (j = 1) and another mode is 
antisymmetric (J = 2) . Taking into account symmetry 
properties, the cross gain coefficient can be expressed as 



312 cos ( 



-4VPiP 2 cos 2 



901P11P2 



a 



dxdy, (8) 



where the integration is made over one of the waveguides, 
and the saturation factor is 

G 4 , = (1 + PrVi 2 + P 2 i/>2 2 ) 2 - 4P 1 P 2 V^V'2 cos 2 0. 



4 



It is evident from ([8]) that the term gi2 cos <fi in this case 
is always non-positive. It turns to zero at Pi = or 
P 2 = 0. 

Since the cross-gain term g\2 cos 4> equally diminishes 
an amount of energy extracted by stimulated emission for 
both modes, the increase of the mode power is favored 
for the mode with higher power. This is rather impor- 
tant conclusion radically differing from the intuitive sug- 
gestion that the mode possessing higher saturated gain is 
dominant at an output of a sufficiently long amplifier. It 
is shown below, that the mode can possess higher modal 
gain throughout the entire amplifier but the part of the 
power carried with this mode in the total power dimin- 
ishes in a course of amplification. 

The idea of using the saturated modal gains [27| stems 
from consideration of the incoherent fields of compet- 
ing modes. Such a situation can be thought of as the 
competition of two signals launched in the fiber amplifier 
from independent sources at the same frequency. The 
wave field intensity is expressed in the incoherent case as 
/ = -PiV-'i + f^I) an d the equations for the two- mode 
evolution are 



dPi 

—r- = guPi, 
az 



dP 2 

—r~ — 922 "2, 

az 



(9) 



gn and 322 are the modal gains of modes 1 and 2, corre- 
spondingly. These equations include so-called cross sat- 
uration of gain associated with the fact that both modes 
saturate gain. In the limit of a weak saturation, equations 
© are reduced to the form closely resembling known dy- 
namic equations for two- mode ring laser ([281]) 



X 
Y 



2X (at 
2Y (a 2 



H2Y- 



O12Y) , 
hiX), 



(10) 



where X, Y are dimensionless intensities of modes 1 and 
2, Oj, and (3j, j — 1,2 are the above-threshold small 
signal gains and self saturation coefficients, respectively. 
812 and 02i are the cross saturation coefficients. If there 
is a difference in the modal gains, then one mode may 
suppress the growth of the other. If the modal gains are 
equal, both the modes lase the same power for the case of 
an inhomogeneously broadened (0i2#2i < P1P2) ring laser 
[28| . For a homogeneously broadened (6*126*21 > /?i/?2) 
ring laser [2^] equations lfTU|) predict that from two modes 
starting at t = the mode with higher power suppress 
completely the second mode. Practically, the operation 
regime of a laser in the last case is random due to the 
spontaneous emission effects [H, [3(| . 

Actually, the analogy between weak-saturation limit of 
© and laser equations (fit)]) is valid only for low-power 
wave fields, i. e. for laser starting from small signal. 



IV. SYSTEM OF TWO ULTRA-THIN PLANAR 
WAVEGUIDES 

The delta- function (S(x)) is a favorite potential well 
in quantum mechanics. In application to optics, delta- 
function well is a mathematical limit of the high-contrast 



thin planar waveguide with arbitrarily small waveguide 
width d and a high refractive index difference An char- 
acterized by a single parameter dAn. In this limit, the 
waveguide supports single mode, the wave field of which 
is almost constant within the core and extends far out- 
side the core. It should be mentioned, that the core 
cladding difference is restricted by the condition An <C 1 
for the scalar model |T]) applicability. Usage of the delta- 
function as a model waveguide allows one to study waveg- 
uide arrays analytically 31] . 

Figure[T]shows the system of two ultra-thin waveguides 
situated at locations x = ±0. For this particular system 
it is convenient to introduce dimensionless transverse co- 
ordinate £ = x/a. Equation §T§ for guided modes of the 
waveguide array reads 



d 2 tp 

~de 



+ {-a 2 + 2k [S (£ - 1) + 5 (£ + 1)]} = 0, (11) 



where k = k 2 adnnAn is the delta-function amplitude, a 
is the wave field attenuation rate outside of the waveg- 
uides, a 2 is eigenvalue characteristic for a given mode. 
As well known [32], such a system supports two guided 
modes shown in Fig. [T] provided the condition 2k > 1 
is fulfilled. The fundamental mode is symmetric 4>S (£) 
(in-phasc mode), and the second mode is antisymmet- 
ric, ipA (0 (out-of-phase mode). The amplitudes of these 
modes attenuate exponentially at £ — > ±00 with corre- 
sponding rates as and a a, while in the space between 
the waveguides ips ~ cosh(as£), and ipA ~ sinh(aA<0- 
The attenuation rates satisfy transcendental equations 



a s = K(l + e~ 2as ) 



OLA 



K 1 



-2a/ 



(12) 



The modal shift of propagation constant 5(3 j = (3j — kriQ 
is expressed as 8(3 j = a 2 / (4L fl ) , L R — kn a 2 /2 is the 
Rayleigh length. Gain in the waveguides can be included 
in the model by addition an imaginary part into k: k = 
k! + in" , where k" = gK 1 / (2kAn). The modal gains in 
the small signal limit read 



2a' s K" 
<^s - —j 

Ca 2a' A K" 



sinh 2a' s + (sin 2a'g) a' s /a'g 



sinh 2a' 



S 



cos 2a'g 



sinh2a^ — (sin2a^) a' A /c 



sinh 2a! \ 



cos 2a" A 



where a" and a" are the real and imaginary parts of 



the corresponding parameters satisfying (|12| . For the 
symmetric system under consideration it is convenient to 
introduce a confinement factor of the mode Tj , as being 
proportional to the overlap of the mode intensity and 
gain in one waveguide. For the ultra thin waveguides the 
confinement factor is equal to the squared amplitude of 
the mode in the waveguide. The modal gain is expressed 
as Gj = 2T 1 gd/a. As far as k' > k" (2k An > g) , the 
expressions for modal gains and confinement factors can 
be simplified. The confinement factors of the symmetric 
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and antisymmetric modes can be expressed as 



[l + exp(-2a' s )] 2 
2 l + 2K / exp(-2a^)' 

K 1 [1 - cxp(-2a^)) 2 
2~ ' 1 - 2«/exp(-2a^)' 



Ts = - 



T A = 



ctj (j = S, A) satisfy equations (fT2j) with k replaced by 
k'. In the limit of weak coupling between waveguides 
(k' 3> 1) the antisymmetric mode has higher confinement 
factor, than the symmetric one: 

Ts« y [l + 2«'exp(-2a' s )r\ 

As coupling between waveguides increased the antisym- 
metric mode tends to transform into a leaky mode, 
(a A — > 0) , and its modal gain diminishes proportional 
to a a, while the symmetric mode gain remains to be of 
finite value. Thus, there is a critical distance between the 
waveguides, at which both gains are equalized. The con- 
finement factors are shown in Fig. [5] as functions of n' . 
The curves intersect at k' ci — 0.900126 .... This value 
corresponds to a ~ 3.9 /jm for the system of two waveg- 
uides with An = 2 ■ 10~ 3 , d = 2 /mi, no = 1.456 and the 
radiation wavelength 2n/k = 1 /im. The Rayleigh length 
Lr in this case is 70 /im. 



V. WAVE FIELD AMPLIFICATION IN THE 
SYSTEM OF TWO ULTRA-THIN WAVEGUIDES 

It is instructive to analyze wave field amplification in 
two ^-function type coupled waveguides taking as a refer- 
ence case the standard equations applicable for descrip- 
tion of two incoherent modes with gain cross saturation 
accounted. For definiteness, the simplest gain saturation 
model is adopted g = go/ (1 + I/Is&t) > and field intensity 
in following is measured in /sat units. 



A. Incoherent wave fields 

For the incoherent fields of two modes in the waveg- 
uide system shown in Fig. [1] the equations for the modal 
powers ((9]) read 



where P a /Pao = (Ps/Pso)\ 7 = T A /T S , Pso and 
Pao are the mode powers at the amplifier entrance. An 
asymptotic (£ — > oo) behavior of the modes depends on a 
value of 7, which is equal to ratio of modal gains. Namely, 
at 7 = 1 both modes grow with equal rates. This case 
is illustrated in Fig. [31 which shows the diagram in di- 
mensionless variables Is = Ps^S and I a = P A T A . It is 
clearly seen that the proportion Pa/ Ps remains constant 
at any distance. If the modes have different overlaps with 
the gain, the mode with higher confinement factor wins 
and its power increases linearly with length, while the 
power of the second mode grows as a fractional power of 
the length. It means that powers of both modes increase 
with no limit. This result does not depend on the initial 
proportion of the powers of the modes, if the amplifica- 
tion is large enough. 



B. Coherent wave fields 

If the wave fields of two competing modes are coherent, 
gain coefficients entering into the system of equations Q 
can be found explicitly from eqns. (0 and © as 



(13) 



g j3 = (1 + P S T S + PaTa) / (CL a ) , 

9as = -4 Ko T A r s /fOTcos 2 4>/ (CL R ) , 

here j = A,S and the denominator contains 

C = (1 + P S T S + P A T A f - 4TaT s PaPs cos 2 



In the specific system under consideration the ratio of 
the modal gains of the two modes is independent from 
the field in the amplifier: 



gss/g 



AA 



r s /r, 



The reason is that both the modes are equally distributed 
over two waveguides and the mode profiles inside the 
waveguides are considered to be identical due to small 
waveguide thickness. The equations © read 



dP s 2k"P s T s 



d( 



C 



dP A 2k"P a Ta 



d( 



C 



(l + P s r s - J P A r A cos20), (14a) 
(l + P A r A -PsTsco82(t>), (14b) 
^*L R 60. (14c) 



dP s 
d( 
dP A 

d( i + p s r s + p a t a ' 



i + P s r s + PaT a 

2k' 'P a Ta 



Direct integration of (|T4j) was performed by a Mathcad 
software solver based on 4 th order Runge-Kutta method. 
The calculations were made for the cases: a) n' = 
0.900126 (T A = T s ), and b) k' = 1.2 {T A = 1.112r s ). 
It was taken for calculations n' ' = knogoda/2 = 0.001. 
This value can be achieved at the small signal gain 
go 0.21 cm -1 for a = 5.2 /im and the other param- 
eters of the construction as at the end of section IIV1 so 
that k' — 1.2. The spectral separation of the modes in 
In P s / Pso + Is (Ps - Pso) + Ts (Pa - Pao) = 2r S K 'C, this case is 6/3 = Q.277^ 1 w 22 cm" 1 . The results of 



where C — z /Lr and 



go (Lad/ a) . This system of 



equations can be easily integrated 
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calculations are presented in Fig. [4] For the equal con- 
finement factors = Tg behavior of curves on diagram 
in fig. 4a showing I a as a function of Is for different 
initial conditions is in contrast with their behavior in the 
case of incoherent modes (Fig. [3|. In the last case, both 
modes increase power so as the proportion of powers is 
constant. In the former case, there is the line Is = Ia, 
which is unstable in Lyapunov's sense: an infinitesimal 
deviation from this trajectory results in further running 
away curves approaching asymptotically either vertical 
or horizontal line in dependence on direction of initial 
deviation. This behavior corresponds to dominance of 
one mode, power of which grows linearly with the length, 
while the power carried by the other mode is stabilized at 
a certain level. This behavior is similar to one observed 
in the Lamb's model [28| for two-mode laser in the case 
when the gain cross-saturation coefficient is greater than 
the self-saturation coefficient. 

If the confinement factors are different (see Fig. 4b), 
the diagram on the whole is nearly the same. The straight 
line Is = I a in the diagram Ia(Is) in fig- 4a transforms 
to a curve (not shown in Fig. 4b), the shape of which 
could be found numerically. This curve is a separatrix di- 
viding (Is, I a) plane into two parts. In the upper part all 
curves approach vertical asymptotes (asymmetric mode 
dominates), while in the lower part all curves approach 
horizontal asymptotes (symmetric mode dominates). 

To get better insight into mechanism leading to de- 
pression one of the modes in two-mode amplification let 
us illustrate the general arguments adduced in section ITiTl 
by results of analysis of the specific construction under 
consideration. Figure [5] illustrates behavior of the terms 
in the right hand side of eqs. ([6|) gjjPj (j — A, S) and 
VPlPsgAS cos 4> with modal gain and cross gain coef- 
ficients defined by (H3J). Pso = 0.27 and P A o = 0.13. 
Small-scale oscillations of these terms associated with 
mode beating are shown in Fig. 5a. It is seen that an- 
tisymmetric mode power increase qaaPa is lower than 
that of symmetric mode despite that its modal gain is 
higher. The cross gain component V 'PaPs9as cos 4> is 
non-positive and antiphase to the power increments asso- 
ciated with modal stimulated emission. The same terms 
averaged over oscillations are shown in Fig. 5b. For con- 
venience of presentation the cross gain-term is multiplied 
by (—1). It is clearly seen in Fig. 5b that decrement in 
emitted power caused by the cross-gain term tends to 
equilibrate the term qaaPa at a sufficiently long ampli- 
fication length. In dimensional variables the Rayleigh 
length at parameters taken is 123 ^m. This phenomenon 
leads to stabilization of the antisymmetric mode power 
on the level reached to that moment. 

The trend to domination of one of modes in simulta- 
neous amplification of two coherent modes can be illus- 
trated by behavior of the mode power ratio Pa/Ps- The 
proportion Pa/Ps calculated for k' = 1.2 is shown in 
Fig. [6] as a function of propagation distance for various 
input proportions at the total input power kept constant 
Pao + -Pso = 0.4. is seen that increasing Pao/Pso re- 



sults in change of amplification regime from dominance 
of symmetric mode to asymmetric one. It is worth to 
note that the sign of derivative at C = of the proportion 
Pa/Ps cannot serve as a criterion for the change of the 
amplification regime. As far as Pao/Pso grows, the curve 
appears, for which growth at small distances changes to 
falling down at longer distances. Figure [6] proves that 
there exists a critical value of Pao/Pso, which separates 
regimes of dominance of symmetric or asymmetric mode. 
Actually, the critical value of Pao / Pso is a function of the 
total power and depends parametrically on confinement 
factors values. 

A series of calculations allows us to find the criti- 
cal fraction of the symmetric (in-phase) mode in the 
launched signal shown in Fig. [7] as a function of the 
total power for confinement parameters values (r^ = 
1.112r,s). At such values of Ta and Ts the small sig- 
nal gain of the antisymmetric mode is greater than of 
the symmetric one. Despite this fact, the in-phase mode 
dominates at the output of an amplifier when its fraction 
in the total launched power is above curve shown in Fig. 
[7] The higher is the total launched power the less is the 
critical fraction of the in-phase mode. It follows from 
this figure that 40% excess of the in-phase mode power 
in the input signal is sufficient to suppress the out-of- 
phase mode power for P > 0.3 The analysis performed is 
strictly valid for the model system under consideration. 
However, the mechanism of weak mode suppression by a 
strong mode is of quite general nature. Thus, it is ex- 
pected that this mechanism works in any fiber amplifiers 
provided the input signal has a sufficiently narrow spec- 
tral width in order that results found for monochromatic 
wave field were applicable for a real signal. 



VI. CONCLUSIONS 

Theory of monochromatic wave field amplification in a 
waveguide array is developed. An approach based on ex- 
pansion of the wave field in terms of guided array modes 
leads to appearance of new terms in a system of ordi- 
nary evolution equation for the mode amplitudes. These 
terms have meaning of cross-modal gain and, as shown, 
completely change behavior of the amplified wave field. 
Analysis of two-mode amplification reveals new unusual 
features in characteristics of two-mode amplifier. Instead 
of unlimited growth of both modes for incoherent fields, 
the effect of weak-mode suppression by a strong one takes 
place. Detailed analysis is made for the amplifier com- 
posed of a pair of ultra-thin waveguides. The critical val- 
ues of waveguides and input signal parameters are found, 
at which the in-phase mode dominate at the amplifier 
output. The conditions for asymptotically stable single- 
mode amplification are found. The model developed is 
applicable for studies on multimode competition in the 
multi-core fiber amplifier. 
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FIG. 1: Schematic of the system of two parallel ultra-thin waveguides. The profiles of the symmetric (solid line) and antisym- 
metric (dashed line) modes correspond to k' = 1.2. 



1.0 n 




FIG. 2: Confinement factors for the symmetric (solid line) and antisymmetric (dashed) modes vs. the coupling strength 
parameter k' . The quantities are dimensionless. 
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FIG. 3: Mode intensities in the waveguide normalized to /gat for different values of initial powers ratio for the incoherent wave 
fields of the modes, Fg = Ta- 
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FIG. 4: Mode intensities in the waveguide normalized to 7s a t for different values of initial powers ratio: a) Ts = Ta; b) 
Fa = 1.112rs- Dashed line in fig. 4(b) is bisector Ia = Is- 
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FIG. 5: Dimensionless terms qas cos (f>\/PAPsLn(so\id line), qaaPaLr (dashed line), and gssPsLn (dotted line) as functions 
of propagation distance normalized to Lr: (a) axial oscillations shown within a short propagation interval; (b) averaged over 
oscillations values, gAS cos 4>\/P~aPsLr (solid line) is taken with minus sign. Pso = 0.27 and Pao = 0.13; Fa = 1.112rs- 
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FIG. 6: Proportion of the modal powers for varied inputs at constant total input power as a function of the amplifier length. 
The quantities are dimensionless, Fa = 1.112r,s. 
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FIG. 7: Critical fraction of the in-phase mode in the input total power P — Pao + Pso, above which this mode dominates. The 
quantities are dimensionless, Fa = 1.112Fs. 



